Transverse energy distributions and J/ip production in Pb+Pb collisions 
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We have analyzed the latest NA50 data on transverse energy distributions and J/ip suppression 
in Pb+Pb collisions. The transverse energy distribution was analysed in the geometric model of 
AA collisions. In the geometric model, fluctuations in the number of NN collisions at fixed impact 
parameter are taken into account. Analysis suggests that in Pb+Pb collisions, individual NN col- 
lisions produces less < Et >, than in other A A collisions. The nucleons are more transparent in 
Pb+Pb collisions. The transverse energy dependence of the J/ip suppression was obtained following 
the model of Blaizot et al, where charmonium suppression is assumed to be 100% effective above a 
threshold density. With fiuctuations in number of NN collisions taken into account, good fit to the 

■ data is obtained, with a single parameter, the threshold density. 

C^l I PACS numbers; 25.75.-q, 25.75.Dw 

In relativistic heavy ion collisions J/ip suppression has been recognized as an important tool to identify the possible 
phase transition to quark-gluon plasma. Because of the large mass of the charm quarks, cc pairs are produced on 
a short time scale. Their tight binding also make them immune to final state interactions. Their evolution probes 
the state of matter in the early stage of the collisions. Matsui and Satz predicted that in presence of quark-gluon 
plasma (QGP), binding of cc pairs into J/ip meson will be hindered, leading to the so called J/ip suppression in heavy 
CnI ion collisions Q. Over the years several groups have measured the J/ip yield in heavy ion collisions (for a review of 
^ the data and the interpretations see (|,^). In brief, experimental data do show suppression. However this could be 

attributed to the conventional nuclear absorption, also present in pA collisions. 
^ID , The latest data obtained by the NA50 collaboration Q on J/-0 production in Pb+Pb collisions at 158 A GeV is the 
' first indication of anomalous mechanism of charmonium suppression, which goes beyond the conventional suppression 
' in nuclear environment. The ratio oi J/ip yield to that of Drell-Yan pairs decreases faster with Et in the most central 
,—1 . collisions than in the less central ones. It has been suggested that the resulting pattern can be understood in a 
C3 ' deconfinement scenario in terms of successive melting of charmonium bound states 

In a recent paper Blaizot et al ||] showed that the data can be understood as an effect of transverse energy 
• fluctuations in central heavy ion collisions. Introducing a factor e = Er/ETib) and assuming that the suppression is 
I ' 100% above a threshold density (a parameter in the model) and smearing the threshold density (at the expense of 
another parameter) best fit to the data was obtained. Capella et al |^ analysed the data in the comover approach. 
(— I ' There also, the comover density has to be modified by the factor e. Introduction of this adhoc factor e can be justified 
• • . in a model based on excited nucleons represented by strings 

At a fixed impact parameter, the transverse energy as well as the number of NN collisions fluctuate. The Fluctuations 
in the number of NN collisions were not taken into account in the calculations of Blaizot et al |^ or in the calculations 

■ of Capella et al In the present paper, we present a calculation following the model of Blaizot et al ||] which includes 
5^ I these fluctuations. As will be shown below, if the fluctuations in number of NN collisions are taken into account very 

good description of the NA50 data can be obtained without smearing the threshold density. The smearing effect is 
generated by the fluctuations. 

Geometric model has been quite successful in explaining the transverse energy production as well as multiplicity 
distributions in A A collisions ||,^. In this model, it is assumed that at impact parameter b, the number n of NN 
collisions is Poisson distributed with average < Ub >. In the Glauber approximation < Ub > is written as,, 

< Ub J d^sTAis)TBis - b) (1) 

where is the inelastic NN cross-section, assumed to be 32 mb. In this model all the nuclear information is 

contained in the nuclear thickness function is Ta_b(s) = ^ dzpA,B{s, z). In the present calculation, we have used the 
following parametric form for pA{f) 

l + exp{—^) 

with a=0.53 fm, tq — I.IA^/^ . The central density is obtained from / pA{r)<fir — A. 
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In the geometric model, the probabihty to obtain Et at impact parameter b in n- number of NN collisions is written 
as, 



L{n+ I) 

where Q^^^{Et) is the n-fold convolution of Et distribution resulting from elementary NN collisions, 

Q^'^^Et) = J dE^...dEJ^g{E^)...g{EJ^) x 5{Et - E^... - EJ^) (4) 

In eq.^ 9{Et) is the normalized Et distribution for NN collisions. Most of the Et distributions in NN collisions 
can be well approximated by the gamma distribution, 



g{ET) = ^e-"^-i?^-i (5) 



with the parameters a and [3. For the gamma distribution, the average and the variance are, 

<ET>NN=l3/a (6a) 
<El> / <Et>^ -l = l/f3 (6b) 

Gamma distribution has an elegant convolution property which greatly facilitate computation, n-fold convolution 
of a gamma distribution is again a gamma distribution, with parameters a' = a and f3' = n(3. Thus, 

Q^''HET)^^^e-'^''-Ef-' (7) 

The final transverse energy distribution is then obtained from eq.^by summing it over n (from 1 to oo) and averaging 
over the impact parameter b. 

Geometric model has been quite successful in explaining the Et distributions in heavy ion collisions with a ^ 2 and 
/3 ^ 2 Q . We have fitted the Et distribution in Pb+Pb collisions varying the parameter a and /?. NA50 collaboration 
did not correct the Et spectra for the efficiency of target identification algorithm, which is lower than unity for Et 
lower than 60 GeV. To obtain the parameters a and /3 we have fitted the purely inclusive part of the Et spectra 
{Et > 60 GeV). Very good fit to the data is obtained with a = 3.46 ± .19 and /3 = 0.379 ± 0.021. The fit is shown 
in fig.l. We note that best fitted values of a and /3 indicate that average Et produced in individual NN collisions is 
order of magnitude smaller than in other AA collisions Q. The variance is also a order of magnitude large. It seems 
that the nucleons are more transparent in Pb+Pb collisions than in (say) 0+Au collisions, leading to less < Et > in 
individual NN collisions. It is evident that Et production mechanism in Pb+Pb collisions is different from other AA 
collisions. Transparency of nucleons in Pb+Pb collisions may be interpreted as an indication of QGP production in 
Pb+Pb collisions. 
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As mentioned in the beginning, we have followed the model of Blaizot at al |g[ to analyze the Et dependence of 
charmonium suppression. Charmonium production cross-section at impact parameter b is written as, 

(fa^'^/d?b - a-^/^ j dhTX^f{s)T^J^'{s - h)S{h, s) (8) 



where T^^^ is the effective nuclear thickness function, 



T^ff{s)^ dzp{s,z)exp{~aabs dz'p{s,z')) (9) 



with aabs as the cross-section for J/^ absorption by nucleons. The exponential factor is the nuclear absorption survival 
probability, the probability for the cc pair to avoid nuclear absorption and form a J/ip- S(h, s) is the anomalous part 
of the suppression. Blaizot et al jsj assumed that J/ip suppression is 100% effective above a threshold density (ric), a 
parameter in the model. Accordingly the anomalous suppression part was written as, 

5(b,s) = eK-np(b,s)) (10) 

where Up is the density of participant nucleons in impact parameter space, 

np(b,s) = Ta{s)[1 ~ e-^""^«(b-^)] + [Ta ^ Tg] (11) 

Recognizing that the endpoint behavior of charmonium suppressions are due to transverse energy fluctuations, 
Blaizot et al modified the density of participant nucleons by a factor e = Et / EtQ)). This modification makes 
sense only when rip is assumed to be proportional to the energy density. Implicitly it was also assumed that the Et 
fluctuations are strongly correlated in different rapidity gaps. The assumption was essential as NA50 collaboration 
measured Et in the 1.1-2.3 pseudorapidity window while the J/4>^s were measured in the rapidity window 2.82 < y < 
3.92 1^. In the geometric model, strong correlation between Et fluctuations in different rapidity window is readily 
obtained. At a fixed impact parameter, fluctuations in Et can be calculated as 

< E^ >AA - <Et >\a ^ <N^ > - <N >^ 1 <E^ >NN - < Et >%n 

<Et>\a " <^>2 <iV> <Et>%j^ ^ ' 

The Et fluctuations has two parts, (i) geometric in nature which remains same irrespective of rapidity window, (ii) 
which depend on fluctuations in the Et distributions in NN collisions. The 2nd part changes with rapidity window 
but its effect is less as it is weighted by the factor 1/ < N >. Fluctuations of Et in different rapidity windows are 
thus strongly correlated. 

We calculate the J/ip production as a function of transverse energy, at an impact parameter b as, 

oc 

da-'/^/dET = Y,Pn{h,ET)P{il^ I EtM (13) 

n=l 

where Pn{h,ET) is the probability to obtain Et in n NN collisions (eq.(||)) and P(V' | ET,\i) is the probability to 
produce a charmonium with transverse energy Et- Piip I ET,h) is given by eq.^, with anomalous suppression part 
modified according to, 

5(b,s)-e(n,--^np(b,s)) (14) 
np/a 

where we have replaced ET^b) by nf3/a appropriate in the geometric model. This modification takes into account the 
fluctuations in number of NN collisions at fixed impact parameter b. 

The Drell-Yan production was calculated similarly, replacing Piip \ i?T,b) in cq.^ by the Drell-Yan production 
cross-section, 

d^a^^/d% = a^^ J dhTAis)TBis - b) (15) 
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c7abs=6-4 mb, nc=3.S fm"2 
(Tabs=4.0 mb, nc=3.44 fm"2 
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FIG. 2. J/?/> survival probability in Pb+Pb collisions as a function of transverse energy. 



In fig. 2, we have compared the theoretical charmonium production cross-section with NA50 experimental data. 
The normalization factor a'^l^ ja^"^ was taken to be 53.5. The solid curve is obtained with (Tafcs=6.4 mb, and = 
3.8 fm~^. Very good description of the data from 40 GeV onward is obtained. It may be noted that if the fluctuations 
in the NN collisions were neglected, equivalent description is obtained with threshold density ric = 3.75/to^, with 
smearing of the Q function at the expense of another parameter. It is evident that in this model, the smearing is done 
by fluctuating NN collisions. Theoretical calculations predict more suppressions below 40 GeV, a feature evident in 
other models also. It is possible to fit the entire Et range, reducing the J/^'-nucleon absorption cross-scetion. The 
dashed line in fig. 2, corresponds to aabs^4: mb and ric — 3.42/to^^. 

To summarize, we have analysed the transverse energy distribution in Pb+Pb collisions as well as the charmonium 
production data obtained by the NA50 collaboration. The transverse energy distribution was analysed in the geometric 
model. It was seen that in order to fit the experimental data, individual NN collisions are required to produce (on 
the average) less E^, compared to other AA collisions. It seems that in Pb+Pb collision, the nucleons become 
transparent. The charmonium production data was analysed following the model of Blaizot et al including the 
effect of fluctuations in number of NN collisions at fixed impact parameter. The experimental data from 40 GeV 
onwards could be very well fitted with (Tabs=6.4 mb and a threshold density of 3.8 fm~^. Neglecting the fiuctuations 
in number of NN collisions, equivalent fit could only be obtained by smearing the O function at the expense of an 
added parameter. Data in the entire Et range could be fitted reducing the J / ^-micleon absorption cross-section to 
o'abs=4 mb. The threshold density is also reduced to 3.44 fm~^. Considering that nucleons become more transparent 
(as suggested by the Et data), such a reduction seems plausible. Melting of charmoniums above a threshold density 
as well as apparent transparency of nucleons (as evident from analysis of Et distribution) strongly suggests that in 
Pb+Pb colhsions, QGP like environment is produced. 
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